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A change of the particle density (by gas pnff, pellets or imparity seeding) dnring 
the plasma discharge in tokamak prodnces a radial cnrrent and implicitly a torqne 
and rotation that can modify the state of conhnement. After ionization the newly 
born ions will evolve toward the periodic neoclassical orbits (trapped or circnlating) 
bnt the hrst part of their excnrsion, which precedes the periodicity, is an effective 
radial cnrrent. It is short, spatially hnite and nniqne for each new ion, bnt mnltiplied 
by the rate of ionization and it can prodnce a snbstantial total radial cnrrent. The 
associated torqne indnces rotation which modify the transport processes. We derive 
the magnitnde of the radial cnrrent indnced by ionization by three methods: the 
analysis of a simple physical pictnre, a nnmerical model and the neoclassical drift- 
kinetic treatment. The results of the three approaches are in agreement and show that 
the current can indeed be substantial. Many well known experimental observations 
can be reconsidered under this perspective. In reactor-grade plasma the conhnement 
can be strongly inhuenced by adequate particle fuelling. 
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I. INTRODUCTION 


Every event of ionization of a neutral particle in tokamak plasma it followed by the 
displacement of the newly born charges (electron and ion) towards the equilibrium orbits. 
They leave the magnetic surface where they have been created, due to the neoclassical drift, 
and evolve towards stationary trajectories. Since the electrons drift much less than the 
ions, the main effect is associated with the newly born ions. Neglecting collisions, the ions 
will settle on circulating or trapped orbits and during the periodic motions they depart 
radially, relative to a certain magnetic surface, alternatively to larger and respectively to 
smaller radius. Since these positive and negative radial deviations relative to the magnetic 
surface compensate, the time average shows no effective radial displacement: the orbit has 
an effective “center” which corresponds to the spatial average of the successive positions 
of the ion (for example: the “center” of a banana, for the trapped ion; we neglect smaller 
neoclassical motions of this “center”, like the toroidal drift). However there is a part in the 
radial excursion of the new ion which remains uncompensated. This is precisely the first 
interval, just after the ionization, when the ion evolves to take the periodic trajectory, and its 
successive positions do not yet average to the “center”. This displacement, from the position 
where ionization takes place, towards the “center” of the periodic trajectory, is an effective 
radial current. At the end of this finite, transitory part, the motion becomes periodic and 
there is no radial current. The radial current of the first part is a source of torque, implicitly 
rotation and in this way has an impact on the quality of the conhnement. It effectively makes 
a connection between the density change (via pellets, gas puff or impurity seeding) and the 
change of the conhnement. We note that there is a considerable experimental evidence 
that the variation of density (during the discharge) produces a change in the quality of the 
conhnement. 

On a fast time scale a radial electric held is generated by the charge separation: the 
electrons are almost tied to the magnetic surface, while the ions will travel with the neo¬ 
classical drift velocity to their “center”, on a distance about half of a banana width. We 
estimate the radial current and obtain an order of magnitude of the rate of the torque. 
Compared with the damping rate of a poloidal rotation by transit time magnetic pumping, 
the ionization-induced rate can be substantially higher. Some improved regimes in JET, as 
“pellet enhanced performance” (PEP)^,— , DIII-Di^,- and conhnement changes observed 
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in many devices^,—,—,— appear to be connected with this effect of density variation. 

The ionization-induced torque has a direction which is hxed by neoclassical orbit’s ge¬ 
ometry and it interacts with any pre-existing rotation which may have been induced by 
Reynolds stress, Stringer mechanism or by external factors (NBI, ICRH). The new torque 
can enhance the pre-existing rotation or can act against it, which makes difficult to predict 
its consequences in all situations. 

We suggest that this process may be a unifying connection between a wide class of 
regimes where it has been noted a correlation between a dynamic change of density (within 
a discharge) and the change of the conhnement. 

According to the preceding explanation, there are two mechanisms that are responsible 
for this connection 

(1) the change of density via ionization (of a pellet, gas puff, impurity seeding or influx 
of neutral atoms from the edge) means that \fe (e: = rjK) fraction of the newly created 
ions are trapped and are moving radially to occupy the positions (the “centers”) which are 
the averages of positions on the trapped (banana) orbits. While after arriving there the 
bounce averaged radial displacement is zero for the trapped ion, the hrst step, when the 
ion moves from the place where it has been created to the “center” of the banana is a net 
radial current, a single and unrepeatable event for every ionization event. The ensemble 
of such events is a radial current that produces a torque which generates poloidal rotation 
and sustains it against magnetic pumping2i; the sheared poloidal rotation is a barrier that 
reduce the turbulence and enhances the conhnement. 

(2) every conversion of a trapped ion into a circulating one (and equally the reversed 
process), is accompanied by a substantial radial drift. This is because the “centers” of the 
two kinds of orbits are different and the change from one type of periodic motion {e.g. 
trapped) to the other type of periodic motion (circulating) goes through an intermediate 
regime, unique and transitory. It consists of the last part of the motion on banana, when 
the periodicity is lost, followed by the hrst part of the motion, until the new periodicity 
is established. Both these parts are unique and transitory and are manifested as a radial 
current which produces a torque. Then any dynamic process which implies a (slow or fast) 
change in the velocity space in the region: trapped/circulating ions, will produce a torque. 
It is interesting to note that this is a mechanism of direct coupling between the toroidal and 
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poloidal rotation. If for some reason a toroidal flow occurs in plasma, a number of trapped 
ions will have increased their parallel velocity and will change from trapped to circulating^^. 
The associated radial current produces a poloidal torque. Conversely, stopping the toroidal 
rotation leads a subset of the ion population to convert from circulating to trapped, which 
again produces a transient radial current and further poloidal rotation. 

In the present work we concentrate on only the hrst of the two processes. 

Several effects can be connected with this ionization-induced torque. Regimes that are 
conhrmed by experiments, like the mentioned PEP, density peaking and/or anomalous den¬ 
sity pinch may have a connection with this torque. The dynamic charge separation that 
occurs when the ions move from the place of ionization to the “center” of the neoclassical 
periodic orbit induces a return current of the background ions. Since the number of new 
ions generated by ionization of a pellet is episodically comparable with the local background 
ion density, the response motion of this latter population takes the aspect of a massive, 
even if short, radial displacement. On a relatively large space interval, on which the radial 
derivative of the rate of ionization [dS/dx in Eq. fl2^ below] keeps the same sign, the dis¬ 
placement of the background ions has a unique direction and is sustained all along the total 
time of ionization. This appears as a density pinch, eventually contributing to the density 
peaking. Since this equally involves the impurity ions, it can provide a new mechanism for 
the impurity accumulation. 

We note that the ionization of impurity atoms leads to much larger radial drifts and 
in consequence larger radial currents and torque. This must be examined in relation with 
impurity (argon) seeding at the edge and with Li pellets in the core. In general any influx of 
neutral atoms in the plasma will be a source of rotation which affects the local conditions, 
including the possibility of being a trigger for the L to H mode transition. 

We give in the next Sections a simple description of the statistical build-up of a radial 
current associated with the ionization. At this level of description the collisions are neglected, 
as are the dispersions of the absolute magnitude of velocities, and of the parallel velocities. 
The intention is to draw attention to the high amplitude of this current. Further, our result 
is conhrmed by the drift-kinetic neoclassical approach, developed in parallel to the classical 
treatments of Rosenbluth and Hinton for the similar cases of the torque induced by alpha 
particles^^ and by neutral beam injection (NBI)^. 
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II. ESTIMATION OF THE RADIAL CURRENT GENERATED AT 
IONIZATION 

A. The contributions to the current 

For an easier discussion we adopt a simple picture, mainly having in mind the pellet 
injection. The source is considered to be limited to a hnite segment [ri,r 2 ] on the radius r, 
placed somewhere between r = 0 (magnetic axis) and r = a (edge). Due to the symmetry 
we take the segment as lying in the equatorial plane. The words “left” or “right” refer to 
this segment, with left being closer to the magnetic axis. The newly born ions will move to 
place themselves on the periodic neoclassical trajectories, banana or circulating. A radial 
current is produced during only the hrst, unique and transitory, part of the trajectory, which 
is about half of the width of the banana. 

We start with a purely geometric example. Consider the motion of a point on a circle 
C with radius R and center at {x = R,y = 0) on the Ox axis. We introduce the angle 
9 {t) between the Ox axis and the radius from the center of the circle (R, 0) to the position 
of the point on the circle {x,y). 9 increases clockwise. The point starts from the origin 
(x = 0,1/ = 0), 9 {t = 0) = 0. The equations are: x{t) = R — Rcos [9 (t)] and y (t) = 
Rsin [9 (t)]. The motion on C is assumed uniform 9 {t) = ut. The average of the positions 
of the point {x (t) ,y (t)) up to the current time t are x{t) = ^ f* dt'x (F) = R — \R^ sin (cat) 
with T (t = 0) = 0 and y{t) = \ dt'Rsin (cat') = ^ [1 — cos (cat)], with y {t = R) = 0. 
Clearly the asymptotic t —)■ cxa average position will be (R, 0), which is the “center” and, 
after a transient phase, the motion is periodic. The speed with which the average position 
moves is dx/dt = (R/ca) t“^ sin (cat) — Rt“^ cos (cat). We note that the //-projection of the 
average, y (t), is always positive and that for small time, t ca“\ the cc-projected average 
is linear in time x (t) ~ (Rca^/3) t. Both properties will also be found for banana orbits. 

The speed of the ions on this transient part of its motion is the neoclassical drift ve¬ 
locity \Di = (1/f^ci) n X 4- (n ■ V) nj ^ (^'^1/2 + t'jjj /R- Here Dei is the ion 

cyclotron frequency, are respectively the perpendicular and the parallel velocity of the 
ion, fx = v^/ (2R) is the magnetic moment, n = B/ |B| is the versor of the magnetic held 
B and R is the radius of curvature of the magnetic line. We introduce the notation vdi , 
constant and positive. The sign of the velocity will be given according to the particular type 
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of ion’s motion and according to the helical orientation of the magnetic held line. The latter 
is given by the direction of the plasma current, which we take in the following anti-parallel 


to the main magnetic held. 

We consider hrst the trapped (t) ions that have, at the moment of ionization, parallel 
velocity in the same direction with the magnetic held (-1- ||) and we note that their banana 
is entirely outside the magnetic surface on which the ionization has taken place. The ion 
radial displacement, on the length d = A*"*" half of the ion banana and only on the time 
interval of this displacement 6t = A^^/vdi, is toward the edge. 
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Figure 1: The banana orbit of an ion, which is fully inside the magnetic surface where it has 
been created. The magnetic surface is represented by the cyan dashed circle. The open black dot 
corresponds to the center of the positions r (t). The open red dot (and the small red vertical line) 
indicate the center of the positions on the coordinate x (t). Just for better visibility we show the 
orbit dilated along the radius by an arbitrary factor (xlO). 

The trapped ions that have in the point of ionization a velocity anti-parallel to the 
magnetic held vector, (— ||) have banana orbit entirely inside the magnetic surface and the 
transitory displacement of the new ion, until the “center” of this banana, is towards smaller 
r radius, d = — < 0 (where A*“ is dehned positive). 

In addition there are circulating (c) new ions, i.e. on untrapped orbits. Those that have 
at the initial point a momentum directed parallel to the magnetic held (-1- ||) have a circular 
orbit which is entirely inside the magnetic surface on which the ionization has taken place. 
The ehective displacement to the virtual center is toward larger radius, d = A'^’*' > 0 (A'’’*' 
is dehned positive). 
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Figure 2: The banana orbit of an ion, which is outside the magnetic surface. The magnetic 


surface is represented by the cyan dashed circle. The open black dot corresponds to the center of 


the positions r (t). The open red dot (and the small blue vertical line) indicate the center of the 
positions on the coordinate x (t). Just for better visibility we show the orbit dilated along the 

radius by an arbitrary factor (xlO). 


The last type consists of ions that are circulating and with their velocity at the initial 
point anti-parallel to the magnetic field (— ||). For them the closed orbit fully includes the 
magnetic surface and it is equivalent to the displacement of the average position of the new 
ion to smaller r, toward the magnetic axis, d = —A^~ < 0 is dehned positive). 

The four contributions to the current through (r, t), from the two pairs (he.(± ||)-trapped, 
respectively (± ||)-circulating) appear to be small and they partly compensate, having op¬ 
posite signs. However, even the small remaining (effective) radial displacement, multiplied 
with the rate of generation of new ions, leads to a signihcant radial current. There is an¬ 
other aspect: the centers of the two bananas discussed above [(+ ||) and (— ||)] are spatially 
separated and we can associate to them the full population of new ions generated in those 
two positions. The rate of ionization, expressed in number of ions per per second has a 
signihcant radial variation, both for pellets and for gas-puff: the rate of generation of new 
ions may differ substantially between two radial positions, even if they have small separation, 
of the order of centimeters. The two banana centers are separated by a distance A*"*" -|- A*“ 
which is in this range and suggests that the radial variation of the rate of ionization is an 
important factor. 

For any point r, the four contributions combine into a single, short-lived, hnite-spatial 
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Figure 3: The banana orbits of two ions. One is inside (red) the other is outside (blue) the 
magnetic surface. The magnetic surface is represented by the cyan dashed circle. The two open 
black dots correspond to the centers of the positions r (t), for each orbit. The two open red and 
black dots (and the small red and blue vertical lines) indicate the centers of the positions on the 
coordinates x (t) for both orbits. For better visibility we multiply the radius coordinate of each 

orbit by an arbitrary factor (xlO). 

size - event of ion charge displacement, i.e. a current. These events occur everywhere within 
the radial segment of ionization and for all the time when there are still neutral atoms to 
be ionized. 


B. Calculation of the flux of electric charge of the new ions into a point (r, f) 


The pellet contains a total number of particles (neutral atoms) Nt and the ionization 
takes place in the toroidal volume between the surfaces ri and r 2 , Vt = 2 'k^R {r\ — rf) . The 
total time for ionization of the particles of the pellet is The rate of generation of ions 























Figure 4: The untrapped (circulating) orbit of an ion, fuly inside the magnetic surface. The 
magnetic surface is represented by the cyan dashed circle. The open black dot corresponds to the 
center of the positions r (t). The open red dot (and the small red vertical line) indicate the eenter 
of the positions on the coordinate x (t). Just for better visibility we also show the orbit dilated 

along the radius by an arbitrary factor (xlO). 




Figure 5: The untrapped (circulating) orbit of an ion, fully outside the magnetic surface. The 
magnetic surface is represented by the cyan dashed circle. The open black dot corresponds to the 
center of the positions r (t). The open red dot (and the small red vertical line) indicate the center 
of the positions on the coordinate x (t). Just for better visibility we show the orbit dilated along 

the radius by an arbitrary factor (xlO). 

per unit volume and per second has the average magnitude 

. ioniz M, 1 

n {r,t) ~ (1) 

^ ^ ' q-ioniz ^ ^ ' 
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Figure 6: The untrapped (circulating) orbits of two ions fully enclosing (blue) respectively fully 

inside (red) the magnetic surface. The magnetic surface is represented by the cyan dashed circle. 
The two open black dots correspond to the centers of the positions r (t), for each orbit. The two 
open red and black dots (and the small red and blue vertical lines) indicate the centers of the 
positions on the coordinates x (t) for both orbits. For better visibility we multiply the radius 
coordinate of each orbit by an arbitrary factor (xlO). 


The rate of ionization has strong spatial and temporal variation and we hud more convenient 
to express it as 


. ioniz .ioniz 

n (r, t) = Hq S (r, t) 


/ ions\ 

\m^s J 


( 2 ) 


where S (r, t) > 0 is by dehnition a nondimensional function representing the space-time 
variation of the ionization rate. The maximum of S' is 1 and its variable is normalized, 

. ioniz 

X = r/a. The constant factor Hq (of order ~ ions/m'^/s) is the physical quantity 
that carries information on the average rate of ionization and is taken from experimental 

. ioniz 

observation. The factors Uq and S are constraint by the condition 



. ioniz 


dVn 


ir,t) = Nt 


( 3 ) 
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Figure 7: The time evolution of the average radial position r (t) for an ions whose orbit is fully 
inside the magnetic surface. The straight horizontal line is the asymptotic average position Vasymp 
(the “center” of the periodic orbit). The dashed vertical line in the asymptotic range is used to 
select the part of the trajectory to determine Vasymp, far from the transient. The open dot at the 
intersection of r = Vasymp and r = r (t) determines the approximative end of the transitory 
regime, Tasvmv ■ The drift velocity is estimated as vdz = where rini is the radial 

^ ^ '^asymp 

position at ionization. 


For example, taking for simplicity a source that is constant in time for the entire duration 
of the ionization [0,r“"’*^], we have 


|•r2 


Un 


S (r) A'K'^Rrdr = Nt 


'ri 


( 4 ) 


The hnite volume Vt is divided into toroidal shells of infinitesimal width dr on the minor 
radius, placed at the position r , with volume Vdr = ^dr = Air'^Rr [dr) where A is the surface 
area. We fix a reference position r G [ri, r 2 ], and calculate the net flux of ions that traverses 
the surface A = Air'^Rr at r. Consider the bunch of ions that are produced in a time 5r, 
filling the elementary shell (denoted D) situated at a distance p from r. Their number is 

. ioniz / \ 

n (r — p,t — j X Vdrdr. The ions from D that are generated at time t — p/vci travel 
with constant velocity voi and arrive in r at time t , traversing the surface A in a time 6t. 
The flux at {r,t), in {ions/m?/s), is 


F (r, t) = n 


p \ 11 .i 

r - p,t - VdrOT—— = n 

VDi J AdT 


- pR- —) 

VDi J 


dr 


( 5 ) 
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Figure 8: The time evolution of the average x-projection of the position, i.e. x{t), for an ions 
whose orbit is fully inside the magnetic surface. The details are the same as for Figure [71 



Figure 9: The time evolution of the average y-projection of the position, i.e. y{t), for an ions 
whose orbit is fully inside the magnetic surface. We note that it remains positive at all times. 


The maximum distance of travel on which a new ion generates a current is from the point of 
ionization until the “center” of the periodic orbit, pmax = A (which is one of A*^). Between 
r — A and r there are many inhnitesimal shells, at distance p' from r (A > p' > 0). The ions 
created in these intermediate cells arrive at time t in r if they are generated at t — p'/voi- 
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Figure 10: The time evolution of the average radial position r (t) for an ions whose orbit is 
outside the magnetic surface. The details are the same as for Figure [71 



Figure 11: The time evolution of the average x-projection of the position, i.e. x (t), for an ions 
whose orbit is outside the magnetic surface. The details of the Figure are the same as for Figure 

m 

Summing these partial contributions Eq.((5]), the origins of which are in the interval [r — A, r] 
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Figure 12: The time evolution of the average y-projection of the position, i.e. y (t), for an ions 

whose orbit is outside the magnetic surface. 



Figure 13: This plot represents the time evolutions of the average positions on r, i.e. r{t), for 
the trapped particles whose orbits are fully inside (lower curve) and respectively outside (the 
upper curve) the magnetic surface. The details are the same as in Figure [71 


14 







Figure 14: The time variation of the average positions on r, i.e. r (t), for the two types of 
trapped orbits are represented here with the purpose to give an idea of the smallness of the 
difference in the average radial displacements. One of the curves (the lower one in Figure flBI) has 
been reversed in order to make easier the comparison. 



Figure 15: The same as Figure [TT] but for the ^-projection. 
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Shells of infinitesimal dr width, sources of ions for the final (r) position 



Figure 16: A simple representation of the fluxes that are traversing the surface placed at the 
radius r, coming from various infinitesimal regions 5r. The sum over these contributions is the 
integral that defines the current density in (r, t) according to the text. The positions ri and r 2 

are the limits of the domain of ionization. 


it results F^ (r, t) 


. / . lomz / , O \ . 

T^{r,t)= n (r — p,t - dp 

Jo V '•^Di 

pA 

. lomz I 

^ no dp 

Jo 


( 6 ) 


= rin 


nr ^ fdS\ , fdS\ p' 
, 1 fdS\ ,2 1 fdS\ 


where the derivatives of S are calculated in (r, t). This flux must be multiplied with the 
fractional number representing how many of the new ions will settle on trapped, respectively 
circulating orbits. We take the approximative values ^/e and respectively l — ^/e. In addition, 
we assume that a fraction of 1/2 new ions have parallel, respectively anti-parallel initial 
velocities. 
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C. The current density in {r,t) 


Taking into account the four type of ion’s orbits, we use Eq.([H]) to estimate the flows of 
new ions coming in, or leaving, the point {r,t). The contributions from neighbor points are 

/■A*+ 


h™ (r,i) = ^v^ 


1 


n 




VdiJ 


dp' 


(7) 


1 . ioniz 

= 




dS\ {A'+y 


dt J VDi 


This is the number of ions that are trapped and had an initial velocity parallel with B. They 
are produced in (^r — p',t — , summed over [r — A*"*", r], and flow toward (r, t) from the 

left (he. their current is positive). 


r!l ir,t) = -ve 


. lomz 

n 


r + p',t 


P 


dp' 

VdiJ 


( 8 ) 


1 .ioniz 

2V^^o 


S{r,t) A*- + 


dr J ^ 2 \dt J VDi 




This is the number of ions that are trapped and had an initial velocity anti-parallel to 
B. They are produced in + p',t — , summed over the interval [r, r -|- A*“] and flow 

towards (r, t) from the right (he. their current is negative). 


r'h"’ (r,() 


1 

2 


(i-vi) 


I 



. ioniz 


n 


0 


^ (r, t) A‘’+ 






(Ac+)2 

VDi 


(9) 


This is the number of ions that are circulating and had an initial velocity parallel with 
B. They are produced in {r — p'A — ^ j) summed over the interval [r — A'^+, r] and flow 
towards (r, t) from the left (he. their current is positive). 


rT">(r,«) = i(i-v7) 


. loniz 

n 


T ^ p'd — — ] dp' 

Vdi. 


( 10 ) 


= 5(1-^/?) 


. loniz 
ric 


2 \ dr J ^ 2 \dt J VDi 


This is the number of ions that are circulating and had an initial velocity anti-parallel to 
B. They are produced in (^x + p',t — j, summed over [r, r -|- A'’ ] and flow towards (r, t) 
from the right (he. their current is negative). 
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The current density resulting from the in fluxes is 


j‘" {X,t) = |e| (r‘i+"> - r',',-*' + r'g" - r'i""') (ii) 

or 

r^{x,t) (12) 

7 nvi') y 

= -\e\n, [S (x, t) + {l - ^) A'^^ - {l - ^) A^] 

+5 (f) [-^ Ah' - (1 - V?) (A-)= - (1 - Vi) (A-)^; 

+5iy (f) [-A Ah' + ^/^(A'A - (1 - Vi) (A-)^ + (1 - Vi) (A'-)^]} 

We now have to calculate the flows {out) that leave the shell of the point (r, t).The groups 
that are leaving (r, t) are: 

Trapped with parallel (he. positive) initial velocity t {+ ||) 

x/_i_||\ 1 . ioniz 

^out (toward large r) (13) 

The second group consists of trapped with anti-parallel (he. negative) initial velocity 

in 1 /— • 

^out (toward small r) (14) 

The third group consists of circulating ions with parallel initial velocity 

(-I l|\ 1 , • ioniz 

Uo),/ (r, f) = - (l — \/e) n {r,t) (toward large r) (15) 

The fourth group consists of circulating ions with anti-parallel initial velocity 

riout {'i'',t) = - [l — y/e) n {r,t) (toward small r) (16) 

Summing the out components after taking into account the signs according to the de¬ 
scription 

j“‘ (r,«) = |e| (r® - r«;"> + rf - ri») = o (i7) 

Adding the flows in and out J (r, t) = J*” (r, t) — (r, t) we obtain 

J {r,t) (18) 

1 . ioniz , r / , , , X / ___ / , XT 

= 5 kl {s (r,«) [vk (A‘+ - A'-) + (1 - Vi) (A'+ - A-)] 

(f) [-A (Ah' + Ah') - (1 - ^kF) ((A-)' + Ah'); 

(w) [A (- (^‘h'+Ah') + (1 - A (- (A-)'+Ah')]} 
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We show in Appendix A that the contributions of the circulating ions is much smaller 
than that of the trapped ions and for the present estimation can be neglected 

1 


, • V L/ 11/1/ >0 ^ ^ t I t ^ 

J(r,f) - |e|no (x, f) [A + - A “] 


(19) 


In the numerical model (next Section) we take a time-independent source, dS/dt = 0, 
and the current density becomes 


J (r) 


1 .ioniz 

- |e| n, 

2 


'-0 


?^S(r,i) (A‘+- A-) - i ^ 


(A‘+)"+(A‘-) 




( 20 ) 


The distance A*^ travelled by the new ion from ionization to the “center” of the periodic 
motion (he. the distance on which there is effective current) will be calculated in the next 
Section by solving the equations of motion of the ion, as initial value problem. For the 
present estimation we adopt neoclassical approximations^ replacing A*^ with the “radius” 
of the banana, 

A ~ 'VjJiThounce PiQ^ ^ 

where ions Voi ^ % and ~ rjve = In regions where the ionization 

rate has strong spatial variation, the second term in Eq. fl20l) is large and we can simplify 
the result as 

(22) 


''0 


2 2 - 1/2 
* 1 ^ 


For an estimation we take a = 1 (m), R = 3.5 (m), Bt = 3.5 (T), R = 1.5 {keV), 
Nt = 3 X 10^^ neutral atoms in the pellet, 4 x 10“^ (s) duration of the complete 

ionization process^,— and the radial extension of the zone of ionization is between ri = 0.4a 
and r 2 = 0.7a. The energy of the new ions is a fraction [t] = 0.75) of the background 
ion energy and the trapping parameter A = h = 1 -|- ecosd is taken A = 0.92. It 

• ioniz 

results Vt ~ 22.8 (m'*), A (r) 138r ( m^) and the average rate of ionization n ~ 3.3 x 

10^^ {ions/irA/s). Adopting for S (r) a simple spatial prohle limited between = 0.475a, 
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Tfo = 0.625a , with max S' = 1 we hnd from the constraint Eq.(jl]) Uq k, 1.2 x 10 
{ions/w?/s). Using values suggested by experiments^^,— A*^ from exact integration, 
and a~^dS/dx ~ 20 (m“^), we obtain from Eq. fl20l) \ J\ ~ 11.52 {A/m?). 
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The result is indeed high. For comparison we consider nrriidve / dt ~ JBt- This means 
dve/dt ~ 23 X 10® {m/s^), or, in every microsecond the poloidal speed would increase with 
more than 4 (km/s). In less than one tenth of a milliseconds vq rises to the range of the 
ion thermal speed Vth,i ~ 0.38 x 10® {km/s). For comparison the transit time magnetic 
pumping decay of the poloidal velocity would contribute with 

= I'"”'= i 0 + 


where I is the mean free path and Pa is the ion-ion collision time^. We hnd ~ 

1.2 X 10^ and taking a poloidal velocity (as observed in some experiments, e.gA) 

ve ~ 10^ (m/s) it is estimated {dve/dt)^^ ~ 1.2 x 10® (m/s^). It results 




> 



(24) 


Of course, this ionization torque acts for short time (few milliseconds) and Eq. fl2^ is an 
overestimation as long as the neutral atoms’ dynamics {e.g. the pellet cloud) and the 
bulk ion’s reaction (return current) are not described in detail. But this result is a strong 
suggestion that the ionization torque is important. 


III. NUMERICAL IMPLEMENTATION 

The numerical simulation of this process has been done on a discrete mesh x 

{ik}k=iNT ~ ~ ^00. The simple physical picture described above has been 

implemented.. For any cell {ri,tk) we calculate the number of ions that are generated, 
Uq S {r,t), using the expression for S {r) = a + (3 {r /a) + 'j {r /a) . Imposing S (r^) = 
S (rfe) = 0 and max S' = 1 at (r^ -f r^,) /2, the coefficients are determined and, as explained 

. ioniz 

above, Eq.(jl]) determines the constant Uq = 1.2 x 10^® {ions/m'^/s). We now assume 
that the energy of the new ions is a fraction (0.75) of the background ion thermal energy 
(at Tj = 1.5A:eU) and that the probability of being trapped is y/e. Finally we assume that 
the probability that the ion has an initial velocity which is parallel to B is 1/2, equal with 
the probability to be anti-parallel. Now we look at the way the ions move. The total 
excursion on r is A*®" (to larger r) and (to smaller r). The displacement is represented 
on the mesh and every cell {i',k') which is traversed by the flux of ions stores 
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this contribution, adding it to a variable that will finally be the current fiowing through it. 
There are several other cells whose new ions traverse this cell {i', k') and all contributions 
are counted and summed. Due to the assumed constancy of there is no decay along 
the path that starts from (rj,ffc) and ends in the cell (r^ + [respectively 

(r, — A*“,f + A*“/n£)j) for the anti-parallel initial velocity]. All cells traversed along this 
path retain the contribution from {ri,tk). 

Profile of the (x=r/a,t)-dependent part S(x,t) of the ionization rate 

0.9^ 

0.8 
0.7 
0.6 
0.5^ 

0.4 
0.3 
0.2 
0.1 ^ 

0 

15 

X 10’'' 

5 

Figure 17: The space-time profile of the ionization source, the function S {r,t). The factor 
no {ions//s) multiplies this function to obtain the effective rate of production of new ions. 



Instead of the neoclassical approximations for A*^ and vdi— we choose to solve the system 
of equations of motion^,-,-,—,- in every cell (r*, Nx-k=i nt- 


dr 

dt 


1 
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dS ^ ^11 
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Radius of the asymptotic position (center) (m) for +V|| 



Figure 18: The (r, t) profile of the final position (the “center”) of the average r (t) for the ions 
whose orbit encloses the magnetic surface. This is calculated by solving in every for 

i = 1, NR, k = 1, NTIME cell of the disctere mesh, the set of equations of motion for a trapped 

ion. 

From the solution we get the exact orbit of an ion born in the cell (r*, NR-k=i nt t’ut 
we still have to operate the separation of the transitory part, the part which represents the 
unique manifestation of a current, from the periodic part of the trajectory, whose average 
does not produce a current. We calculate for each trajectory the time-dependent average 
position r {t) = j r (F) dt' and leave the integration sufficiently long such that the asymp¬ 
totic quasi-static average position Vasymp = f(t^ oo) to be clearly identihed. This position 
^asymp (the “center” of the banana) is retained and the quantity A* is obtained as the dif¬ 
ference between Vasymp and the initial position, which is the point of the ionization, Vini- 
We still need the estimation of the effective time that is necessary for the ion to reach this 
asymptotic position. The hrst intersection between the asymptotic line r = Vasymp and the 
evolution line r (t) takes place at a moment tasympi which is retained as the representative 
time of ion’s travel to the center. This procedure is admittedly approximative but we have 
tried several reasonably alternative methods and the present one seems the best. 

Once we know Vasymp — 'Tini = A and tasymp we hnd voi = ^/tasymp and all data for 
the displacements of the new ions originating from {ri,tk) on the mesh are now available. 
In every cell the system is solved for both parallel and anti-parallel initial velocities and 


22 




FLYP, distance from ionization to the center (m) for +Y| 





time (s) 


Figure 19: The variable FLYP, difference between the final position (the “center” of the 
banana) and the initial radial position (where ionization occurs) for an ions whose orbit encloses 
the magnetic surface. This is the amount of radial displacement of the ion on which there is 

current. 


we calculate the distances the times L^gymp of excursion and the drift velocities 
This is shown in Figure 1 for ions with parallel, respectively anti-parallel 
initial velocity. The asymptotic position Vasymp is obtained by averaging a set of late values 
of r [t), for t beyond the vertical dashed lines. This ensures a good precision of identihcation 
of the “center”. The straight line r = Vasymp intersects the evolution line f (t) is a point 
marked by a open dot, the time tasymp- 

We have adopted a prohle S (r) which is constant in time. Figure 2. Its support is [va, h] 
inside the interval [ri,r 2 ]. The empty radial regions on both sides are necessary because 
the excursions of lengths of the ions generated at the ends of the support must be 
recorded in these regions. The result (the current J (r, t)) depends on time since the process 
that starts at t = 0 (no-ionization) rises slowly by accumulating current contributions, 
before saturation. More interesting is the spatial prohle conhrming that the sign of dS/dr 
is decisive and that the total torque, i.e. integrated over the plasma volume, is zero, as 
expected from conservation of angular momentum and from the fact that no ion is lost 
from plasma in our picture. The current is plotted in Figure 3. In Figure 4 we show the 
mesh-cells that contribute to the current calculated in a reference cell, chosen arbitrarily in 
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Time to reach the periodic orbit (s) for +V|| 



Figure 20: The {r,t) profile of the calculated time for an ion to reach the average position (the 

“center”) starting from the place of ionization. 


Drift velocity of a new ion (m/s) for +V|| 



Figure 21: The (r, t) profile of the calculated drift velocity v^i, as explained in Figure [71 

(r,/) = {ix = 350,itime = 150). The dots aligned on the a straight line (since vm =const) 
at left represent cells from where the new ions with positive (parallel) initial velocity arrive 
in r at t. The straight line at right represents ions with anti-parallel initial velocity. 
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Radius of the asymptotic position (center) (m) for -v,. 



Figure 22: Same as Figure [THl but for ions whose orbits, calculated in every mesh-cell, are inside 

the magnetic surface. 


FLYM, distance of displacement of a new ion (m) for -v 



Figure 23: The variable FLYM, similar with FLYP from Figure [T9l but for the “smaller” 

banana. 
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Time to reach periodicity (s), for-v 



Figure 24: The time to reach the “center”, for the “smaller” banana. 


Drift velocity of a new ion (m/s) for -V|| 



Figure 25: The drift velocity for the smaller banana. 


IV. DRIFT-KINETIC CALCULATION OF THE IONIZATION-INDUCED 
CURRENT DENSITY 


Similar problems are treated in the Ref.— (for alpha particles) and the Ref.— (for NBI). 
For the present case the drift-kinetic equation for the new ions, is 


dl 

dt 


+ (niih -b VDi) ■ V/ 


5 


ioniz 


(25) 
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The mesh-ceils that contribute to the current lor 


The mesh-cells that contribute to the current, for + 



X list of mesh-cells for FLY, -v. 


X list of mesh-cells lor FLY. -t 


Figure 26: At left: the mesh-cells that are traversed by the ions generated in the point 
(ir = 250, itime = 25) that evolve towards smaller radii (their banana orbits are fully inside the 
magnetic surface. At right: the mesh-cells that are traversed by the ions that evolve towards 
larger radii (their banana orbit ecloses the magnetic surface). 


where the drift velocity of the guiding centre is 'Voi = “f'nn x V the neoclassical 

notations will be used: ^ = v\\/v = (1 — A/h)^'^^, A = hv']_/v‘^, h = 1 + ecosO. The limits 
of the trapped particle region in the variable A are 1 — -^<A<1-|--^. The velocity space 
variables are v, A and a (=sign of uy). 

Since in this simple treatment we neglect the effect of collisions the neoclassical small 
parameter is the ratio of the banana half-width to the minor radius, S = A^^/a 1. In usual 

neoclassical perturbative solution of the drift kinetic equation the zero order distribution 
function is the Maxwellian. In the present case, the perturbative expansion of / , the 
solution of the drift-kinetic equation for the new ions, must contain a term which is directly 
related to the source and, since this is determined by external factors, it cannot be ordered 
as powers of S. The hrst term is formally of order —1, /_i. 


/ — /-I + /o + /i + 


(26) 


The lowest order 

ii ■ V/_i = 0 (27) 

shows that /_i is constant along the magnetic lines, or /_i = /_i {ip). This result is con¬ 
nected with an assumption about the distribution of ionization processes in space: they have 
a rate which is constant on a magnetic surface. 
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The cells that contributes to a REFERENCE cell 
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Indices of X mesh-cells 


Figure 27: The mesh-cells that contribute to the current density that is “measured” in a 
reference cell (r, f), with indices of discretization ir = 350 and itime = 150. From the cells that 
are marked by a dot the ions are moving to larger radii (for the dots that are at the left of the 
axis of symmetry of the figure) and respectively to smaller radii (for the dots that are at the right 
of the axis of symmetry). Accordingly the contributions must have been generated by ionization 
in a position r — p and respectively r + p, and at a time t — p/v^i, for them to reach the reference 
cell at (r, f). The real mesh is too detailed to be shown and only two patches are shown, for 

illustration. 


The zeroth order equation is 

«lin . V/„ = -v„ . - ^ + ■S”"" (28) 

As in any multiple space-time scale analysis we average at this level (0) to obtain a solution 
on the level (—1). We apply the operator of bounce averaging to eliminate the function /q. 
This gives the equation 

df-l -^oniz - =—df_i , , 

^ = S - v„ . (29) 

The operator of bounce averaging is A = ^ ~ T I bounce time 

is T = f = f VII/(qR) • limits of integrations for untrapped ions are [—tt , tt] and 

for trapped ions the integral is dehned 




Number of contributions in a single cell, (lX.ITIME) 



0 0.4 


x-coordinate 

Figure 28: The {r,t) profile of the number of contributions that are registered in the cells of the 
mesh, coming from neighbor cells, as described in Figure [271 


the current density (A/m*^) 



Figure 29: The (r, t) profile of the current density J 


where —Oq and + 6^0 are the turning points of the banana. The radial projection of the 
guiding centre drift velocity can be written 

v„.V^. = /„,(n.V)(^)=/„,B.v4(^) (31) 

where I = ■ V(p = RBt = I ("0) is a function of only the magnetic surface variable (V’) 

and h-V 6 ' = 1 / {qR). At this point we assume that the new ion has reached the asymptotic 
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the current density (A/m^) calcualted using the analytic formula 



time (s) 


space(m) 


Figure 30: The first analytic formula that obtains the current density. 


the current density (A/m ) from analytic2 formula 



time (s) 


space(m) 


Figure 31: The second analytic formula for the current density. 


periodic motion on the banana. Then the radial displacements average to zero 


_ 1 r+do 

(vd. ■ VV.) = - 5^ / 

^ J-ea 


dd 




-»o qB.de \a. 


(32) 


Eq. fl2^ becomes 


d f—l -^oniz 


(33) 


and indeed /_i appears as a direct result of the “external” source. The source of new ions 
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Figure 32: 
calculation. 


Comparison Numeric (continuous line) - Analytic (dashed), Analytic2 (.) 



The comparison between the two analytic formulas and the result of the numerical 
The section is made at time that is half the total duration of the ionization process. 


of velocity uq, with direction uo and trapping parameter Aq ia^ 

6{v - Vo) 


b = n 


6a,aoQ (t) 5 (A - Ao) 


TTVn 


(34) 


The Eq. fl5^ simply describes the accumulation of new ions with (uo,cro,Ao) on the surface 
^|J. The motion of these ions toward the banana trajectories and the periodic motion that 
follows must be found at higher orders. Returning to Eq. fl28p we express /o in terms of /_i, 
using ([3l]) 

= (35) 

with the solution 

/ Du \ fl-f ^ (qIA 

+ g{7p,e,t) (36) 




dijj 


where a constant of integration of the operator Vy is introduced. We make few remarks. 
First, the time dependence of /o inherited from /_i will be essential for the presence in the 
theory of the hrst, transitory and unique, part of the trajectory. Further, the hrst term can 
be approximated, using for circular geometry 

df_i (-0) B U|| df-i (r) U|| d 


Be Sr’ 


-I 




:/-i(^) 


(37) 


d'ljj BeVt^i dr ^eddr' 

and we see that the difference between the distribution function /o and that of the previ¬ 
ous level (/-i) consists of a radial shift of the space argument, of the order of the poloidal 
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Larmor radius v\\/VLgc.i ~ pg. This is the same relationship as between the hrst order neo¬ 
classical distribution function relative to the Maxwellian equilibrium distributiort^i. In the 
particular case of trapped particles, the correction needs also to reflect the approxima¬ 
tive relation between the thermal speed and the parallel velocity of ions22, and we have 
v\\l^eci ~ {vth/^eci) = pgs^!"^. Finally we note that g is constant on the magnetic 

lines, i.e. on surfaces, h-= 0. A distribution function for bananas can never be constant 
on the magnetic surfaces because the trajectory stops somewhere. Therefore g must only be 
added to (136|) if we consider circulating particles. For our purpose it is not retained. 


The next step is the equation for the first order fi, which is derived from the equation 
written at zero-order 


^ + ti||n • V/i + VDi ■ V/o = 0 


(38) 


This involves the variation of the hrst order correction function fi along the magnetic lines 

ni|V||/i + .... 

VDi ■ V/o = VD* ■ + VD* • (39) 

oO ofi) 


where x V ( ^) • 


and similarly 




^ - ^11 RR ^ ( ^11 


Returning to the initial expression 


VDi ■ V/o = / 


nil 


d / ^11 ^ ^/o d { v\\\ dfo 


qR id%jj \^ciJ 99 99 \^ciJ 9^! 


The bounce average is 


I'+do 


(VD* ■ V/o) = -I 


d9 n|| 


T J-g^ vii/ (qR) qR 
Here we replace /o with its expression in terms of /_i; 




A (3L] ^ - ^ (3L] ^ 

9-ip Ifici/ 99 99 yffci/ 9fi 


Qr 


9'ip 


(40) 


(41) 


(42) 


(43) 


(44) 
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(45) 


■ V/o) 

— (+—( (^) 
T J_g^ \^^|J \ilcij 09 y^lcij O'lp Oil \VLci) Kytci) d9d'il) 

__a / ^ a/_i (V>) aV-i ] 

09 v^ci/ 0%l) v^ci/ 0^) 09 V^ci/ V^ci/ Oip'^ ) 


The first and third terms cancel. In addition we know from fl271) that /_i is constant on the 
magnetic snrfaces. i.e. the second term is zero. It remains 



(46) 


We can now calculate the radial current, using the distribution functions in orders — 1, 
0, 1. This is obtained from the radial projection of the drift velocity, Eq. flTI]) 



(47) 


The current is projected on the radial direction and the result is averaged over the magnetic 
surface, with the operator (^4) = w~^ f Ard9/B0, w = J rd9 /Bq, 


il 


VV) /) = - |e| I d\< t-iin ■ Vf'j (48) 



(j . = \e 


An integration by parts over 9 has been done. Two terms are absent: (1) the order —1 
distribution function /_i does not contribute due to (HT} : and (2) the order 0 does not 
contribute, due to fl4^ . The first order to have a contribution to this current (averaged over 
surface) is /i. From the equation fl55]l we take the term 


n||n ■ V/i = —— - WDi ■ V/o 


(49) 


and Eq. fHHl) becomes 



(50) 


The surface average operator applied on the second term in the bracket vanishes 



(51) 


This is shown by a calculation analogous to that of Eq. fl5^ , using fl4U]) and fITT]) followed by 
the substitution of 0441) . The current is 



(52) 
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We use the zero-order function /o from Eq. 




and take the time derivative, 
2 


d-^v 




aV-i \ 

d'lpdt I 


(53) 


where, in (IM|) we keep the isotropic velocity space integration d^v = Anv'^dv and introduce 
the factorization ( 12 ]) 


df_i -^oniz -ioniz 1 

= b = nQ b [r,t) -—-0 [V — Vo) 


For circular surfaces 


dt 


/■\ \ \^T 

(jr) ~ - |e| 




<t\, I ^ 


aV-i 

drdt 


(54) 


(56) 


We note again the presence of the poloidal gyroradius, corrected for trapped particles (n|| = 
{r/R)^^^Vth), 

Bl / nil ^ ^ 


Bl VO 


O, 


6 Cl 


2 2 2 -1 
= Pe^ = PiQ ^ 


(56) 


This length corresponds to the radial excursion of the new ion, as dehned in our 

previous approach, Eq. fl2ll) 

After replacing the expression of df^i/dt Eq. fISTD we have 

Ur) = - |e| 


(57) 


q cf d ■ ioniz 1 

dvf,ge-^n„ S(r,t)j^S{v-Vo) 


(58) 


. ioniz dS (r,t) , 9 V 
= - leirin —7^- \Pe^) = 


e n, 


. ioniz dS (r, t) 


{A‘R 


dr “ Sr 

It is understood that pe and further = piqe~^l‘^ are calculated at the velocity Uq. As 
before we replace the A-integration with multiplication with fraction of trapped particles 


. ioniz dS (r, f) 9 9 1/9 

Jr = -7|e|% ^ ' 


(59) 


The constant 7 is a purely neoclassical constant and is calculated, for more general 
conditions, iii^^. It includes the exact integration over the trapping parameter A, which is 
contained in ny at vq hxed. The result is 7 ~ 0.38. 

We note that the analytic structure of our result Eq. (j2^ and of its rederivation in the 
neoclassical drift-kinetic theory, Eq. (l5^ are the same as the expression obtained in the 
treatment of Rosenbluth and Hinton for the current induced by the a particles^^. The 
coefficient in our approximate treatment Eq. fl2^ is 7 ~ 1/2. Figure 5 represents the space 
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dependence of the current J (r, to) obtained from the analysis of the physical picture, Eq. fl20|) . 
from the numerical model and respectively from Eq. fl22l) which is also the result of the drift- 
kinetic approach. The time to is chosen at half the total time interval, to avoid the transient 
after the onset of ionization. 

We note that the conclusion of the mentioned paper, that the rotation induced by the 
creation of alpha particles is insignihcant is a consequence of the very small nuclear reaction 
rate. The equivalent parameter, in the present problem, is the rate of generation of new 
ions, which is three orders of magnitude higher in the case of pellets. 

The two treatments (the simple arguments related with the fluxes of ions and, respectively, 
the drift-kinetic equation) lead to the same result but there is an apparent difference between 
them. In the hrst treatment the separation of the trajectory of a new ion in a transitory 
part, where effective current exist, and a periodic part with no effective radial current is 
the key element that identihes the source of current, torque, rotation. In the drift-kinetic 
approach this separation is not obvious. One would expect an “initial value problem” where 
the distribution function would result as integral over the history of the ion’s motion. This is 
not visible. The Heaviside function of the source is not helpful either: it marks the beginning 
of the ionization process but after that every moment of time is a source of new ions and this 
is not represented. We understand however that the separation is implicitly done through 
the velocity space integration. The current is dehned as j ~ |e| {vDi f (fvf {r,v)). The 
late phase of the ion orbit is periodic and the integral mixes to zero the two-way travels on 
banana, leaving only the hrst, non-periodic, part. 


V. DISCUSSION AND CONCLUSION 

The gas puff, the pellets, the impurity seeding and in general any inhow of neutrals 
into plasma produce a substantial radial current and implicitly a torque that can be higher 
than the magnetic pumping damping. It can be shown that it can also be higher than the 
turbulent Reynolds stress and the Stringer mechanism. We have derived a simple analytical 
expression which is conhrmed by numerical simulation. Furthermore, we have re-derived 
it within the neoclassical drift-kinetic approach. All three methods have close quantitative 
results, as shown in Figure 5. We mention, qualitatively, few possible consequences. 

PEP regimes^^ seem to be connected with the ionization-induced rotation that improves 
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the local confinement by creating effective barriers through the sheared poloidal ffow^. The 
duration of the PEP and the density peaking are compatible with ionization-induced rota¬ 
tion. 

The ionization-induced radial current leads to density peaking, in at least three different 
ways. If the gradient of the rate of ionization of a pellet is negative {dS/dr < 0, higher 
ionization rate close to the plasma center, as for pellets launched from high-field side) then 
the current Eq. fl2^ is directed toward the edge. The bulk ions must move toward the 
magnetic axis to compensate this current. Schematically, we consider a new ion that moves 
a distance A toward the edge and then “stops” (actually it moves periodically on banana). 
An ion of the background must move in opposite direction the same distance A. But in that 
moment another new ion is created at a distance A from this position, closer to the center 
and starts moving toward the edge. Then the background ion must continue its displacement 
toward the center to compensate this new current. While any new ion move a distance A 
then stops, the background ion must continue to move to compensate the small currents. 
Quantitatively the two fluxes are balanced but ions from edge can travel very far toward the 
center. Impurity accumulation in the center can also be produced in this process. 

Second, the rotation produced at ionization is necessarily sheared, i.e. vg = vg (r), for 
two reasons. In the regions of positive and respectively negative radial derivative of the rate 
of ionization {dS/dr ^ 0) the rotation has opposite direction (Figures 3 and 5). In addition, 
the background ions must have a local rotation that is opposite to that of the new ions. The 
rate of extraction of the free energy from density gradients is reduced and the turbulence 
will have shorter radial correlation length. The rate of transport decreases and the peaking 
of the density in the center is enhanced by the smaller density diffusion. 

Third, the shear of the poloidal velocity is actually vorticity uj = dvg/dr and when this 
occurs the Ertel’s theorem ^ _ g imposes a redistribution of density. 

A change of the density at the edge, by impurity seeding!^ or by other strong ionization 
event, must now also be regarded as an electric process, due to the charge separation and 
the radial current of the new ions. It implies that very fast plasma responses should be 
expected^^. This may explain observed fast propagation of perturbations, sometimes called 
“non-local”. Fast increase of the radial electric field is able to determine, as a neoclassical 
effect, the reversal of the toroidal rotation2i,i22^i2S. 
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In conclusion, we have presented arguments that the neoclassical displacements of the 
new ions generated at ionization (of gas puff, pellet, impurity seeding) produce a radial 
current that can be substantial. The current is generated from the first part, transitory, 
unique for any ionization event, of the trajectory: between the ionization and the moment 
where the new ion reaches the stationary periodic motion, trapped or circulating. The 
torque resulting from ionization can be substantial and it can generate internal transport 
barriers. Our perspective on some particular regimes may need reconsideration: Pellet 
Enhanced Performance, density peaking, density pinch, regimes with density higher than the 
Greenwald limit, fast propagation of edge effects, influence of the density on the transition 
to if-mode, connection between density and rotation, reversal of toroidal rotation, etc. A 
more detailed investigation of the ionization-induced rotation is requested, for each of these 
cases. 


A. APPENDIX. COMPARISON OF THE CURRENTS CARRIED BY 
TRAPPED AND RESPECTIVELY CIRCULATING IONS 


The equation of the closed orbit (poloidal projection of the orbit of a circulating ion) is 


{x + arof + y^ = r^ 


(A.l) 


a circle of radius tq that is displaced from the magnetic axis with the amount xq = avQ 
where a = {B/Bt) VDi/v\\ =const.“. We have a -C 1 (since v^i <C uy ~ Vth,i) for circulating 
ions. Therefore the radial displacement of the “center” of the orbit of a new ion that is 
circulating, relative to the center of the magnetic surface where it is created, is small. The 
displacement Xq of circulating ions that are created closer to the edge (at higher tq) are 
larger than that created closer to the magnetic axis. For two ions born in the same point on 
the equatorial plane and with velocities parallel respectively anti-parallel to B the difference 
between these displacements is linear in r. 
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This is indeed very small, due to factor. The closed orbit for parallel initial velocity, 
(-I-), is fully contained inside the magnetic surface, which means that the displacement 
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is positive. The closed orbit for the anti-parallel initial velocity, (—), fnlly encloses the 
magnetic snrface, which means that the displacement Xq ^ is negative. In absolnte valne 
is greater than Xq~\ Note that we here nse the geometric “center of the closed orbit” 
and not the asymptotic valne of the average x (f —)■ cxd). The latter are closer the main axis 
of symmetry, for both sign of the initial velocities. 

The cnrrent carried by the new ions for the short time nntil they access the stationary 


d+) 


Xq Xg 


(-) 


= n'^2a^r 


periodic motion relies on the difference between the displacements 
is directed towards the main axis of the torns. The ratio between the width of the banana 
and the displacement of the center for a circulating particle originating from the same point 
is approximately^^ 




= 4 


X, 



(A.3) 


The second term in Eq. (j2n|l is in general greater than the Erst. We then compare the 
current from circulating ions with only this first term. We hnd that the circulating ions’ 
contribution is smaller than this term, which justihes their neglect adopted in the main text. 
We have to compare the charge displacements, including the densities of trapped (n*) and 
circulating (n'’) particles. We expand, taking as reference position the point r where the two 
ions (with parallel and anti-parallel vq) are born, — A^”) = ^ [A(+) -|- A^”)]. In the 
right hand side, and in all expressions where the difference between A^ is not involved, we 
can approximate A = A^ (r) ^ piq (r) /2. 

^t|A(+) _ A(-)| y^|:ln(gr-V2) 




X, 
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X, 


(-) 


2n‘^a^r 


(A.4) 


We use de/dr ^1/R and (p^g) = 2/S.\/e obtaining the ratio 
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R {pi(l) = r — 
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(A.6) 


(A.6) 


The ratio of the densities of trapped and circulating particles is 


nf 




(A.7) 


and employing Eq. flA.3p 
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Then 
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(A.8) 


In the region of interest [ri,r 2 ] it is sufficient to approximate q ~ 1 + /? {r/af with /? a 
constant. Retaining the dominant term, ^ In ~ we have the estimation 
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t |A(+) - A(-)| 




S+) 


Xq Xq 


(-) 




> 1 


(A.9) 


This justifies the neglect of the current from the circulating ions in Eqs. fllQI) and 


B. APPENDIX. THE MOMENTUM OF THE PLASMA ROTATION 
INDUCED BY IONIZATION 

As mentioned in Section m in the present case the total torque is zero and there is no 
problem of conservation of the angular momentum. For more general S {x, t) the source of 
angular momentum and energy related to the “spontaneous rotation” requires a discussion, 
which we here attempt in general terms. 

Assume a slab-like geometry with the plasma immersed in a static magnetic held B =i?e^. 
From external sources it is applied a transversal electric held E = held. There is a 
motion of the plasma in the direction Gy with the speed vg = E x B/R^, which apparently 
violates the conservation of the momentum along the y direction: before applying E there 
was no momentum along y but after that we hnd plasma moving along y with all its particles. 

There is, of course, no “spontaneous generation of momentum”. The momentum that 
would ensure the conservation and which seems to be missing is actually taken over by the 
external helds B and E. They are acting as an intermediate medium transferring the momen¬ 
tum of the guiding center rriiVE to the external structure (in general coils and condensers) 
that maintains these fields. 

We consider charges of density p (r, t) and currents of density j (r, t) in a limited volume 
V bounded by the surface S. From the Maxwell equation one derives the local balance 

V-R = pE + j xB + £o^(ExB) (B.l) 

with S = ^ where the components of the order-two tensor are2^ 

= (B.2) 
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and respectively 


.. = —BiB. - —-5i.B^ (B.3) 

^ /io ' /io2 ^ ^ 

For a point-like particle of charge |e| with trajectory tq (t), 

p(r,t) = |e| [r - ro (t)] {C/m^) (B.4) 

j (r, t) = |e| V(5^ [r - fq (t)] [A/m^) 


Eq. fIB.lD is integrated over the volnme V 


S n da = + 


^0 


m 


dv IE X Bj 


(B.5) 


'V 


where one notes that the quantity ^qE x B is a volume density of momentum and the last 
term 


= So / dv (Ex B) 


(B.6) 


'V 


represents the amount of momentum of the electromagnetic field inside the volume V. This 
underlies the role of the tensor S: the quantity iS-n is the flux of momentum through the 
element of area da of the surface S. The total mechanical force acting on the particles and 
currents can be expressed as the time derivative of the mechanical momentum inside V 


Then the equation becomes 


jp(e) _|_ p()n) _ ^ pmech 

dt 


[ S-iida = + 

dt' 


(B.7) 


(B.8) 


Anything that changes inside V (held, motion) must be balanced by a reaction from the 
world exterior to V. Now we imagine that S is very large, enclosing the toroidal coils, etc. 
such that the helds have vanished on the boundary. Then 


__'pmech _ pern 

dt dt 


(B.9) 


which expresses in the most clear way the idea of this discussion: any modihcation of the 
mechanical momentum of the charged particle must be balanced by an opposite modihca¬ 
tion of the electromagnetic held. For externally applied (B, E) helds, the momentum (and 
angular momentum) is sustained by a reaction against the sources of the helds, coils and 
condensers. 
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When the electric held is generated by charge separation inside plasma, the momentum 
and the energy of the guiding centers will ncessarly involve the momentum and energy 
inside plasma, besides those external to it. After ionization, the new ions take energy by 
interacting with the background plasma and it is with this energy that they move to settle on 
hnal periodic orbits. Their motion produces the layer of unbalanced ion charge at the edge of 
the ionization region and the resulting electric {E^) held is fully dependent on the energy of 
the new ions. Then the momentum and energy of the plasma moving with ve = E^xB/i?^ 
have their origin in the energy that the new ions could get from the background plasma. The 
rotation v^; is a backreaction, like an inertia. The plasma responds to by polarization 
drift of the background ions. Moving on a distance 5x^ = E^ in the held E^ the 

variation in energy is 5W = |e| x E^ = 6 -^hich is the variation in 

energy of the rotation of the plasma. We can see that the plasma rotation ve is signihcant 
only if the background plasma feeds the new ions with sufficient energy for them to reach 
trapping orbits with substantial i.e. such that the charge separation produces a large 
Eh 

C. APPENDIX. PLASMA RESPONSE TO THE CHARGE SEPARATION 
INDUCED BY THE DISPLACEMENT OF THE NEW IONS 

After ionization the new ions move to take their neoclassical periodic orbit. Between 
the point of creation and the “center” of the periodic motion they carry a transitory, short, 
hnite current. As in the main text we consider the ionization to take place in a volume 
limited between the radii ri (left side, closer to the center) r 2 (right side, closer to the 
last closed magnetic surface). The plasma is considered homogeneous and the ionization 
generates ions that move to the right a distance A* with velocity vei while the new electrons 
can be considered imobile. Then most of the volume between ri and r 2 is neutral but at 
the right end ~ r 2 of the ionization interval it results a layer (denoted I) of positive charge, 
of width A*. This is the source of electric held, resulted from ionization, E\ directed from 
r 2 towards ri. The charge layer I and its held E^ are built up on a time scale 5t = A^/vdi 
in which dE^/dt > 0. The background plasma responds by modifying the Larmor gyration 
orbit from the usual circle to a cycloid (actually the new orbit is a prolate trochoid). The 
deformation of the gyration directly indicates the expected ve = ExB/R^ motion and 
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also the asymmetry of the charge distribution along the new orbit. The asymmetry creates 
a new layer (denoted L) of positive (ion) charge, at the left end, ~ ri and an electric held 
opposite to . This held is sufficient to almost cancel E^ inside plasma, leaving in 
the interior a small i?*"* = E^ — E^, directed to the left, like E^. The asymmetry of the 
modihed Larmor orbit is a manifestation of the polarization drift indnced by the variation in 
time of the electric held E^ (implicitely E^'^'^), a displacement of the backgronnd ions in the 
direction to which points E^ (he. to the left). Since the ionization continnes to accnmnlate 
new ions in the layer /, hence dE^/dt > 0 , the ion’s drift of polarization hlls the layer 
L at the left end, whose electric held E^ (from ri toward r 2 ) continues to quasi-compensate 
E^ inside plasma. 

Qnalitatively, this picture conforms to the concept of return current, which is nniversally 
invoked as the plasma response to any mechanism that is able to prodnce rotation of only 
some component of the density: NBI, ICRH, alpha particle, etc., to which we add: ionization. 
In the following we examine the density of charge and respectively the current density arising 
from ionization, deformation of the Larmor orbit and hnally the polarization drift. The 
ionization-indnced charge separation and current are regarded as “external” factors since 
they are reqnested by the geometry of the held. 

We hrst inclnde a justihcation of the neglect of the volnme-charge accnmnlation that can 
be associated with the strong vorticity. 


A. Charge and current related to the vorticity 

At the edge of the tokamak in the iL-mode regime there is a layer of strong poloidal 
rotation, with radial extension of abont a banana width calculated for the poloidal magnetic 
held. The variation of the velocity magnitnde is very fast in this layer, or, eqnivalently, the 
layer is a concentration of vorticity a; = V xvorci;~ Taking as usnal v = ve = 
we have cn = V x v = or 

a; = (C.l) 

directed along the magnetic held line. The Laplacian of the electric potential A0 is the 
electric charge density and we have the well known sitnation that a vorticity is formally 
eqnivalent to a density of electric charge. If this charge is qnantitatively important, it must 
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be taken it into account together with the currents 


^ + V-J = 0 (C.2) 

where the “charge” = —p^/sQ is the vorticity 

pV = -souB (C.3) 

We can estimate the magnitude of the charge density . If the poloidal velocity has 
a spatial variation from = 0 at the edge of the rotation layer and reaches amplitude of 
~ 10 {km/s) on a radial extension of 10“^ (m) then u ~ lO”^ and this means 

^ 3 X 10“^ {C/m^) (C.4) 


If the formation of this vorticity layer takes place on an interval controlled by the drift 
of the ions then 


5t ~ 


6r 

VDi 


10 ^ (m) 
30 (m/s) 


3 X 10"^ (s) 


(C.5) 


and the time variation of the charge is 

dp^ p^ 3 X 10“^(C/m^) 

dt 5t 3 X 10“"^ (s) 

This must be compared with V ■ J'^. Taking the value estimated in the text ~ 10 {A/m?) 
and a spatial variation on the same radial extension 6r ~ 10“^ (m) we have ||V ■ J'^|| ~ 
{A/m?). This is much higher than the time derivative of the vorticity-charge, so we 
can neglect this latter component of the physical picture. We must remember however that 
~ 10 {A/m?) is obtained for pellets while in other cases (e.p. neutrals penetrating from 
the edge) can be orders of magnitude smaller. In addition the interaction between ions and 
neutrals in this region is complex-. Then we have to check the possibility to neglect the 
vorticity-charge. 



B. The charge accumulation and the current induced by ionization 

The rate of increase of the density of charge p^ by influx of the new ions in the region of 
unbalanced charge at the right end ~ r 2 of the segment of ionization (the charge layer I) is 
dp^/dt ~ |e| hioniz [0 {f — ^ 2 ) 0 {'f '2 + A* — r)] {C/m^/s), where 0 is the Heaviside function. 
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The width of the layer I is A* and the time scale to £11 with newly born ions is 5t 
Using the Gauss law 


dE^ 

dx 




(C.7) 


we take the time derivation and integrate over x (actually x = r and we use x instead of r 
to underline the ID geometry assumed here), 


dE^ 



^1 ^ioniz 



(C.8) 


According to the source of ionization (neutrals penetrating from the edge, pellets, etc.) 
\\dE^/dt\\ can vary over an interval of three orders of magnitude 10®...10^2 The 

static magnitude of E^ can be obtained from the surface charge density = p^A*, as 
= a^/sQ. 

The current density induced by ionization is calculated in the main text 


Ji (r) Si -1 |e| Ar" f (C,9) 

However this calculation was adapted to a particular class of cases and it must be reconsid¬ 
ered for other cases. We just mention the order of magnitude ~ 10 [A/m^]. 


C. The charge accumulation and the current produced by the ion’s drift of 
polarization 


An ion in a constant magnetic field B = SzB performs the Larmor gyration in the 
transversal plane xOy, on a circle of radius pi with frequency VLd = \e\B/mi. When a 
constant electric field E = is added, the circle is deformed into a curve of the cycloid 
type. Integrating the equation 


dD 


e\E^ lei H dr ^ 

-H“- - X 

TTli TTli at 


(C.IO) 


we obtain 


.X 1 \e\E Vyo 1 / 1 |e|E 

Gi- rrii Dro ° Gw V Gw m,- 




1 \e\E 1 
t+TT- 


1 lelE 


Gw \ Gw TTli 


+ cos (Gwt) -h ^ sin (Gwt) 
sin (Gwt) + ^ cos (Gwt) 


(C.ll) 

(C.12) 
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Gyration and Throchoid orbit of a charged particle 



coordinate X (m) xtO^^ 

Figure 33: The deformation of the pure gyration orbit (red) into a trochoid (blue) under 
the effect of an electric field. A very high value of the electric field (350 kV/m) was used in 
order to make more visible the deformation. The plasma center is at right. 

This curve (see Figure [33]) is a prolate trochoid. The electric E x B motion is along the 
negative y axis. We adopt initial conditions (xq = —pi, yo = 0, = 0, Vyo = Vth,i) that 

are identical for the static 7 ^ 0 as well as for pure gyration (i? = 0). In this way we can 
see how the trochoid is different of the Larmor circle. The orbit has, broadly, two unequal 
lobes. This asymmetry makes that the “center” of the positions of the particle to be shifted 
relative to the one of the pure Larmor gyration. For ions the shift is in the direction of the 
electric held E^. The ions are now more frequently present to the left of the symmetry axis 
of the previously symmetric (circle) orbit. There is an effective concentration of ions at the 
end of the interval on r to which points the electric held produced by the “external” 
(generated by ionization). For electrons there is a shift to the opposite direction but much 
smaller and will be neglected. 

Therefore for a given E^ there is an excess of ion charge at the left (~ ri) end of the 
ionization domain. A new layer (called L) of positive charge is generated at the left end, 
opposite to the ionization-induced layer I. The width 5x^ is the amount of deformation 
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relative to the pure Larmor gyration orbit, i.e. the distance between the center of the 
prolate trochoid and the center of the Larmor circle, when both trajectories start from the 
same initial conditions, but with E ^ 0 respectively E = 0. To hnd it, we calculate the 
time evolution of the averages 


x{t) 



X {t') dt’ 


yif) 



y {t') dt' 


(C.13) 


For large t 


X (t ^ oo) 


I \e\E Vyo _ 1 

rrii ilci ^ciB 


(C.14) 


and 


_ /, \ "^xO . 1 

y{t^oo)=yo- — + - 

U Uf‘'l ^ 


\ ) ^B 


(C.15) 


The center of the new orbit has a shift 6x^ = x{t^ cx)) = -^p-^E. The electric field that 
occurs in the above equation is the internal held i?*"*, i.e. the ionization-induced held E^ 
from which we substract the held generated by the new layer L, = E^ — E^. Using the 
shift Sx^ = the surface charge density in the layer L is 


a = e n 


bg_ 


VtriB 


-E 


int 


c 




The electric held produced by the deformation of the Larmor gyration is 


a i 
E^ = — = — 
^0 ^0 


e n 


bg_ 


l^int _ 


^ciB 

from where we hnd = E^ — E^ = E^ — 


-E' 


int 


or 


^int _ 


E^ 


(C.16) 


(C.17) 


(C.18) 


1 c^/v\ 

If the “external”, ionization-induced, electric held E^ continues to increase, there is in¬ 
crease in time of the deformation of the trochoid 


—X (t ^ oo) 
dt ^ ^ 


1 

QciB dt 


(C.19) 


which is precisely the polarization drift of the ions 
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(pol) 

Di 


1 dU®®®* 

QciB dt 


(C.20) 


i.e. 


the drift of polarization simply consists of the time variation of the deformation v 
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Di — 


ix (f ^ oo)- 
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The velocity of the polarization drift of the background ions (of density 


,,(po 0 _ 

^Di — 


dE^ 


£o dE^ 
ilciBl + c^/v\ dt |e|n^^ dt 


(C.21) 


is in general much smaller than the hrst order drift ve and than the neoclassical drift 
Voi- The estimated magnitude varies between ~ 0.02...10 (m/s) according to Uioniz is 
determined by slow gas input or pellets. We would be tempted to expect a slower response of 
the background ions. The build-up of the charge layer induced by ionization is 6t = A^/v^i 
while the build up of the charge layer induced by polarization drift is 
However these two time scales are identical. Using Eqs. flC.191IC.20p and flC.lSp 


\ dt 


-'Di 


-1 


(C.22) 


and inserting Eq. fIC.Sp and 5t = A^/voi we find 


= 5t 


(C.23) 


Using again Eq. flC.8H we obtain 
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(pol) 

Di 


1 1 
QciB 1 + c^/v\ 
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^0 
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This can be translated in the language of currents. By definition 


(C.24) 


j(p°0 = |e| 


(C.25) 


and 


= \e\ A^hioniz (C.26) 

and Eq. flC.24l) shows that the polarization current jfp°h jg equal and opposite to the “ex¬ 
ternally” imposed current . 


Jl ^ j(pol) 


(C.27) 


The plasma response j(P°h jg the return current, involving the background ions. 
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